Coupling between phonons and intrinsic Josephson oscillations in cuprate 

superconductors 
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The recently reported subgap structures in the current-voltage characteristic of intrinsic Josephson 
junctions in the high-T c superconductors Tl2Ba2Ca2Cu30io+5 and Bi2Sr2CaCu208+i are explained 
by the coupling between c-axis phonons and Josephson oscillations. A model is developed where 
c-axis lattice vibrations between adjacent superconducting multilayers are excited by the Josephson 
oscillations in a resistive junction. The voltages of the lowest structures correspond well to the fre- 
quencies of longitudinal c-axis phonons with large oscillator strength in the two materials, providing 
a new measurement technique for this quantity. 



The transport properties of highly anisotropic cuprate 
superconductors in c-direction can well be described 
by a stack of Josephson junctions formed by non- 
superconducting material between adjacent supercon- 
ducting copper oxide multi-layers jjj . Recently the obser- 
vation of subgap structures in the current- voltage (/-V)- 
characteristic of intrinsic Josephson junctions in the high- 
Tc superconductors Tl 2 Ba2Ca2Cu 3 Oio +5 (TBCCO) and 
Bi 2 Sr2CaCu20 8+<5 (BSCCO) has been reported [g-g. 
Each individual branch of the 7-V-curve shows a struc- 
ture which can be traced back to the /- ^-characteristic 
of one single Josephson junction in the resistive state. 
These structures seem to be an intrinsic effect, as they 
have been observed both in step edge junctions (TBCCO) 
and mesa-type stacks (BSCCO) of different sizes. The 
characteristic voltages are completely independent of 
temperature and 

external magnetic fields, which rules out any relation 
to the superconducting gap, vortex flow or the thermal 
excitation of quasiparticles. 

It was shown that the pattern of one junction can 
be described phenomenologically by a resistively shunted 
junction (RSJ) model by assuming ad hoc a special struc- 
ture for the current-voltage characteristic of the quasipar- 
ticles [||. It was argued that such a structure might re- 
sult from peaks in the quasiparticle density of states due 
to Andreev reflection between normal and superconduct- 
ing regions. Several alternative approaches, including 
the modulation of the tunneling distance due to Raman- 
active phonons, have been mentioned in S, but up to 
now all suggestions failed to explain the main features of 
the effect. In this paper we want to discuss a different 
mechanism involving phonons by assuming that the lo- 
cal electric field oscillations produced by the Josephson 
effect in a single junction excite infrared active c-axis 
phonons. In the following we will present a simple model 
where we couple the non-linear current-phase relation of 
one junction to a local oscillator. The analytical and 
numerical solution of this model provides a very good 
quantitative explanation of the experimental data. It is 



shown that the peaks in the subgap structure of the dc 
current-voltage characteristic correspond to zeros of the 
dielectric function of the barrier material, i.e. to longi- 
tudinal optical phonons. 

In the minimal version of the RSJ model the total cur- 
rent (density) 



j c shi7 + <jqE + D , 



(1) 



through one of the junctions is the sum of the Josephson 
current, the (ohmic) quasiparticle current / qp = a^E and 
the displacement current density D, where the (gauge 
invariant) phase difference 7 is related to the electric field 
E in the barrier of thickness b by 



hj = 2eEb 



(2) 



Further time dependencies of a microscopic model are 
thereby neglected for simplicity ||||. 

The displacement current D contains the polarization 
P of the barrier medium, D = t^E + P = e^eE. In the 
case of high frequency Josephson oscillations in the range 
of phonon frequencies it is important to keep the full 
frequency dependence of e(u>) or to treat the polarization 
P as an additional dynamical variable. 

Here we assume that the polarization P = nqz is due 
to a c-axis displacement z of ions with charge q, mass 
M and density n in the insulating barrier between the 
copper oxide (multi)-layers. For the motion of the ions 
we assume a simple oscillator 



z + uj z + rz = —E, 



(3) 



which is driven by the electric field E in the barrier. In 
this model the contribution of the oscillator to the dielec- 
tric function is given by 



£ p hM = 1 + 



nq 



1 



£qM uJq — uj 2 + iroj 



(4) 



1 



It is useful to introduce normalized quantities: First 
equ. (|j) and (|^) are divided by the critical current den- 
sity j c , then a characteristic time variable r = tui c with 
lo c = (2e/h)V c is introduced, where V c — RI C — bj c /ao is 
the characteristic voltage determined experimentally by 
the voltage on the resistive branch at the critical current. 
Then we obtain 



(^0(71) - ^2(71)) cos 70 = v 2 /3 cS (v)ji, 



j = sin 7 + 7 + (3 c j - 
A7 = p + n 2 p 



PP ■ 



(5) 
(6) 



with the normalized polarization p = Puj c /j c , frequency 
£1 = ujo/uj c , friction p — r/to c and the McCumber para- 



meter (3 C — RCuj c — uj 2 /ujj, where 



J J 



2ebj c /(he ) 



is the square of the Josephson plasma frequency and 
C = eoF/b is the capacitance of the barrier. The coupling 
constant A = ujf OD /oj 2 = S(uj 2 /lu 2 ) can be expressed by 
the ionic plasma frequency to 2 on = nq 2 /(Meo) or the os- 
cillator strength S of the lattice vibration. In the present 
case the ratio between the phonon frequency ojq and the 
Josephson plasma frequency tuj is large; therefore even 
a small oscillator strength S = u>q) can lead to a 

sizable coupling. 

The equations (||) and (0) for the phase and the lattice 
displacement can be solved numerically to yield 7(f) as 
a function of the current density j through the junction. 
Taking a time average of the phase the dc current-voltage 
characteristic is obtained. 

From the numerical results it turns out that both the 
phase 7(f) and the polarization p(t) oscillate primarily 
with one frequency which is in agreement with general 
expectations for the RS J-model with large f3 c Js) . There- 
fore it is satisfied to neglect the higher harmonics in the 
ansatz 



7 = 7o + vt + 71 sinwi, 

p = Po + Pl COs(ldt + if) . 



(7) 
(8) 



Here v is the time averaged phase velocity, which corre- 
sponds to the dc voltage v =< V > /V c in the resistive 
state. 

The different Fourier components of the equations of 
motion are obtained by the Bessel function expansion 

00 

sin 7(f) = Jn(li) sin(7o + vt + nwt) . (9) 

n=— 00 

The Josephson current contributes to the dc current only 
if v + nui = 0. 

For the fundamental harmonic we obtain u> = v. Us- 
ing this ansatz in the differential equations (||) and (||) 
we obtain a set of equations for the amplitudes of the 
different Fourier components. As equ. (^|) is linear in the 
polarization p, we can eliminate the polarization p and 
get an equation of motion for the phase oscillation alone: 



3 = v(l + ttct(u)7i), 



(10) 



Ji(7i) SU170 = ~-va{v)-f 2 



with 



/3 e s{v) = Pc - 
a(v) = 1 + 



\(v 2 - n 2 ) 



{v 2 - n 2 ) 2 + v 2 P 2 ' 

\v 2 p 
(v 2 - U 2 ) 2 + v 2 p 2 



(11) 
(12) 

(13) 
(14) 



The functions (3 e g(v) and ct{v) are related to the real 
and imaginary part of the phonon dielectric function 
e P h(w) = ei + ie 2 via 



Asff(u) = Pc^l, 
a(v) = l + v/3 c e 2 . 



(15) 
(16) 



Thereby the poles (zeros) of e(u>) are given by the eigen- 
frequencies wxo (^lo) of transversal (longitudinal) op- 
tical phonons for wavevectors k = 0. In this way the 
formalism can be easily extended to an arbitrary number 
of phonon branches and more complicated lattice dynam- 
ical models. 

As the numerical solution of equ. (0) and (|J) shows 
that 71 < 0.1, the equations ©-(jl|) can be linearized 
in 71 and an analytical formula for the /-^-characteristic 
can be obtained: 



j(v) = v + Aj(v) — v + 



1 



2» a 2 + (vM 2 

1 t ( L_ 

v -lm 

2v 2 \(3 eS + ii 



(17) 
(18) 



Higher orders in 71 have been calculated analytically, but 
have only a negligible influence on the I-V characteristic. 

With the help of these relations some special points 
of the /-^-characteristic j(v) near the subgap structures 
can be identified: 

For small voltages v -C ^to = fi = wto / w c one has 
f3efi(v) ~ Pc and o~(v) « 1 and the model reduces to the 
conventional RSJ-model. For this it is well known || 
that there is a voltage jump in the /-^-characteristic at 

At the resonance v — £Ito of the phononic oscilla- 
tor both the effective quasiparticle conductivity a(v) and 
the effective McCumber parameter [3 e g are strongly en- 
hanced and the I-V characteristic in equ. ( |l8| ) reduces to 
the purely ohmic term. This corresponds to a pure quasi- 
particle tunneling current across the junction, while the 
supercurrent and the displacement current are compen- 
sating each other. 

In contrast to this, equ. (|l^) indicates a resonance in 
j(v) near the zeros of /3 e ff, i.e. the eigenfrequencies f^LO 
of longitudinal optical phonons. 

The difference Aw := f^LO — ^max between JIlo and 
the actual maximum v max of j{v) can be estimated as 
Av < 2%, which is almost independent of the choice of 



2 



parameters. Physically this point is connected with an 
oscillation of the external electric field E and the polar- 
ization p with vanishing displacement current density D. 
With the help of equ. ( |lq ) also an analytical formula for 
the intensity 



Aj(^lo) 



1 



2^LOCT(riLo) 



(19) 



of the subgap structures can be derived. Thereby it 
turns out that a small damping parameter p, correspond- 
ing to a weak coupling of inter-layer ions in neighbour- 
ing contacts or equivalently the small dispersion of the 
phonons in c-direction, is crucial for the existence of a 
hysteretic region; no maximum of j(v) can be found for 
P > Pcrit = ^lo — Oto- Also note that the intensity 
Ai max = j c ■ Aj max in ordinary units is proportional to 
the critical current j c , as it has been reported in experi- 
ments ||. 

In addition to this, an analytical expression for the dif- 
ferential resistivity dv/dj can be derived, which is plotted 
in fig. [p. Note that there exists a region of negative differ- 
ential resistivity for voltages v slightly larger than Hlo, 
which cannot be reached in a current-biased experiment 
with a continously increasing (decreasing) bias current. 
The divergence of the differential resistivity at the subgap 
structures might be connected with experimental indica- 
tions for a significantly enhanced noise production near 
the maxima of the /-^-characteristic. 
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FIG. 1. Analytical solution for the differential resistivity 
for realistic TBCCO-parameters given in fig. 3. 

The numerical solution of the differential equations 
(|J) and (|J) with a Runge-Kutta-algorithm is in excellent 
agreement with the analytical results given here, but it is 
difficult to investigate the region of the /-^-characteristic 
with negative differential resistivity with this method. 

From general results for small capacity j3 c || it is to be 
expected that the maximal current j max and the voltage 
Umax of the peak is reduced in the presence of thermal 
noise. Consequently, the vanishing of the derivative j'(v) 



near Olo is rarely observed in experiments or numerical 
calculations. 

In order to reproduce the experimental data in full 
detail, the shape of the quasiparticle characteristic / qp 
has to be modified. Instead of the ohmic term / qp = ero7 
in equ. (Eh an exponential behaviour 



Iqpd) = cxp 



7-1 

Vb 



for BSCCO and a semiconductor-like dependence 



^qp(7) = 



21 



l + ex P (^) 



(20) 



(21) 



for TBCCO have been used with an appropriate value 
of Vb successfully. The main consequence is a non-linear 
background Iq P (v) in the /-^-characteristic, while the 
modifications of the quasiparticle conductivity <j(v) are 
negligible. The shape and the position of the subgap 
structures are almost independent of the choice of / qp (7) 
and the parameter Vb can be determined very well from 
parts of the /-^-characteristic away from the subgap 
structures. 

For an appropriate choice of parameters these ana- 
lytical and numerical results are in excellent agreement 
with the experimental data, as can be seen in figures (|^) 
and (^). All experimental investigations of phonons in 
TBCCO § and BSCCO [§ agree that infrared active 
c-axis phonons are observed in the frequency range of 
the subgap structures. Theoretical calculations j|] show 
that the dispersion in c-direction is small, which corre- 
sponds to a small ratio p/(Olo — Oto) in our model. 
But considerable discrepancies in the published phonon 
data do not allow us to use the values given in the lit- 
erature directly as input parameters. This is also due 
to the fact that usually only transversal frequencies are 
determined in optical experiments or model calculations, 
while the subgap structures should rather be compared 
with longitudinal branches. Also the precise values of the 
McCumber parameter (3 C and the coupling constant A in 
TBCCO and BSCCO are unknown. Therefore the exper- 
imental /-^-characteristics have been fitted with the un- 
constrained parameters /? c ,f2j, A^, pi, which turns out to 
be quite a sensitive method for the determination of these 
quantities, as an optimal fit seems to be possible only in 
a very restricted parameter range. As shown above the 
position of the peaks is given by the zeros of the phonon 
dielectric constant, which is only dependent on the TO- 
frequencies Oj and the ratios Xi//3. The absolute values 
of Xi and (3 C can be used to tune the relative strength of 
the peaks, while the values of pi are important for the 
overall curvature of the structure. 

The best fits for TBCCO and BSCCO are given in 
figures (||) and (||). The prediction for the eigenfrequen- 
cies t'Lo = ^Lo(k — 0)/(27r) of the longitudinal optical 
phonons with lowest frequencies in both materials are in 
ordinary units: ^lo,i=3.65 THz and i>lo.2— 4.70 THz for 



3 



TBCCO with a characteristic voltage V c = 27.1mV and 
i/ LO ,i=2.96THz, i/ LO , 2 =3.90THz and j^lo. 3=5-71 THz 
for BSCCO with V c = 21.8mV. These values are com- 
patible with the results of most optical experiments WM . 
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FIG. 2. Experimental, analytical and numerical 7-V-cha- 
racteristic of TBCCO for C = 375, v b = 0.29, fii = 0.25, 
Ai = 8, pi = 0.03, Q 2 = 0.34, A 2 = 3.5, p 2 = 0.015. 
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FIG. 3. Experimental, analytical and numerical 7-V-cha- 
racteristic of BSCCO for C = 800, v b = 0.337, Oi = 0.267, 
Ai = 10, pi = 0.008, 0,2 = 0.345, A 2 = 19, p 2 = 0.0045, 
tt 3 = 0.468, A 3 = 39, p 3 = 0.02. 



To summarize, the recently discovered subgap struc- 
tures in the /-^-characteristic in the intrinsic Josephson- 
effect in high-T c -superconductors are explained by a cou- 
pling of Josephson oscillations to optical c-axis phonons 
within a modified RSJ-model. This is - to the know- 
ledge of the authors - the first time that a detailed the- 
oretical explanation of this effect has been given. Apart 
from perfect reproduction of the experimental data for 
both BSCCO and TBCCO new insights in the physi- 
cal interpretation of the peak structures are obtained. 
Above all, the peak position is identified with the eigen- 
frequency of the longitudinal optical phonon, which pro- 
vides a natural explanation for the crucial, experimental 



result of the complete independence of the peak posi- 
tion on temperature, magnetic field and the geometry of 
the probe. In contrast to this, the position of the sub- 
gap structures is expected to depend on pressure and the 
isotopes contained in the material. This fact also sug- 
gests a new approach for a direct measurement of this 
quantity, which is usually hard to determine in optical 
experiments. It also turned out that the width of the 
structure is closely connected with the LO-TO-splitting 
of the phonon branch. Also the observed proportionality 
of the intensity of the resonance to the critical current 
can be understood within the model presented above. 

In principle, similar structures in the current-voltage 
characteristic can be expected near the zeros of the corre- 
sponding dielectric constant, if other kinds of excitations 
with a dipole moment are coupled to Josephson oscilla- 
tions between the layers in a similar way. Extensions of 
this work to a microscopic theory of the Josephson ef- 
fect within the tunneling Hamiltonian formalism in the 
presence of phonon bands Cl(k z ) are currently being in- 
vestigated. This might help to clarify the reason for the 
stability of a local oscillation despite of the finite coupling 
of oscillating ions in neighbouring junctions. 
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